[xl, vl, z, y, x] = fun_element_l(xO, vO, @quadrupoleBl, ... 

gB, q, dt, m, pi, nl, p2, n2, ye, nt) 



Output: 



xl(l), xl(2), xl(3) end co-ordinates of particle, on exit plane of quadrupole 

vl(l),vl(2),vl(3) end velocity of particle, on exit plane of quadrupole 

x(l),x(2),... x value of endpoints of arcs along particle path 

y(l),y(2),... y value of endpoints of arcs along particle path 

z(l),z(2),... z value of endpoints of arcs along particle path 



Input: 



xO = [x0x,x0y, xOz] start co-ordinates of particle 

can be some distance before quadrupole 

vO = [vx, vy, vz] start velocity of particle 

gB gradient of quadrupole (T/m) 

q charge of particle (C) 

m relativistic mass of particle (kg) 

dt time interval of each arc (s) 

nt maximum number of arcs to use in this function 

ye =[yex,yey,yez] reference direction of quadrupole 

pl=[plx,ply,plz] point on entrance plane of quadrupole 

nl=[nlx,nly,nlz] normal to entrance plane of quadrupole 

p2=[p2x,p2y,p2z] point on exit plane of quadrupole 

n2=[n2x,n2y,n2z] normal to exit plane of quadrupole 



Require arc function files: 

fun_element_l .m 

fun_drift.m 

fun_arc_2a.m 

fun_arc_3b.m 

quadrupoleBl.m 

parallel.m 

mag.m 



